Creating and probing macroscoping entanglement with light 
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We describe a scheme showing signatures of macroscopic optomechanical entanglement generated 
by radiation pressure in a cavity system with a massive movable mirror. The system we consider 
reveals genuine multipartite entanglement. We highlight the way the entanglement involving the 
inaccessible massive object is unravelled, in our scheme, by means of field-field quantum correlations. 
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Entanglement is currently at the heart of physical in- 
vestigation not just because of its critical role in setting 
the mark between the classical and quantum world but 
also because of its exploitability in many quantum in- 
formation tasks 0. So far, theoretical and experimental 
endeavors have been directed towards the demonstration 
of entanglement between microscopic systems, mainly for 
the purposes of information processing and manipula- 
tion 0. Nevertheless, the possibility of observing non- 
classical correlations in systems of macroscopic objects 
and in situations close to be classical is very appealing 
and efforts have been made along this direction [3- The 
interest has been also extended to micro- and nanome- 
chanical oscillators, which have been shown to be highly 
controllable and represent natural candidates for quan- 
tum limited measurements, quantum state engineering 
and for testing decoherence theories 0, . This inspired 
us to study an optomechanical device as a macroscopic 
system with readily achievable non-classicality. 

In spite of these exciting progresses, it is in general 
still difficult to infer the quantum properties of a macro- 
scopic object. For the case at hand, the properties of the 
mirror are not directly accessible and one has to design 
strategies to infer its dynamics [3, Isl] . Here we discuss a 
protocol to unravel the quantum correlations established 
in a cavity with a moving mirror. Our scheme uses an 
ancillary cavity field interacting with the optomechani- 
cal device. We treat the whole field-mirror-field system 
as intrinsically tripartite and investigate the behavior of 
entanglement between the subparties. We see signatures 
of mirror-field entanglement in the quantum correlations 
between the two cavity fields, which can be quantified 
through a simple reconstruction algorithm. Our study, 
together with weak assumptions concerning the underly- 
ing model, paves the way to probe information about a 
system by getting entangled with it. We provide an as- 
sessment about entanglement inference in a situation of 
current experimental interest 0] where one of the parties 
is not accessible. 



The model. - The setup we consider consists of two op- 
tical cavities labelled a and h, each in a Fabry-Perot con- 
figuration, sharing a movable mirror. The input mirrors 
of the cavities are assumed to be fixed and each cavity is 
driven by an external field of frequency ojij , input power 
Pj {j — a,b) and coupling strength Ej. The system is 
sketched in Fig. ^ The field of cavity j, locked at the 
frequency ujj ~ uij , is described by the annihilation (cre- 
ation) operator j (j^). In terms of field quadratures, j = 
{xj +iyj)/^/2. The mirror is modelled as a single bosonic 
mode with frequency Um and mass /i. It undergoes quan- 
tum Brownian motion due to its contact to a bath at 
temperature T given by background modes. For stan- 
dard Ohmic noise characterized by a coupling stength 
7m , this leads to a non-Markovian correlation function of 
the associated noise operator ^ of the form {^{t)^{t')) = 
(7m/wm) / u;e-*"(*-*')[l -hcoth(/3a;/2)]dw/27r (with (3 = 
h/ksT and ks the Boltzmann constant) 0- As shown, 
under realistic conditions of weak coupling of the mir- 
ror to the environment, a Markovian description can be 
gained. In this setting the mirror motion is damped at a 
rate 7™. 

In a frame rotating at the frequency of the lasers, the 
energy of the system is written as 

H = ^if^^)+n ^ [(Ao, - Go^m + iE0 - j) 

3=a,b 

(1) 

Here, p, q are the mirror dimensionless quadrature opera- 
tors, Goj = (-l)^^"Goj with G^, = {u:,/lj){h/{ixuj^)f'^ 
the optomechanical coupling rate between the mirror 
and the j-th cavity (length ij), Agj — luj — luij and 
\Ej\ = {2KjPj /{hujij)y^'^ with Kj the j-th cavity decay 
rate. 

In order to study the evolution of the system we refer to 
the Heisenberg picture. The intrinsically open dynamics 
at hand is well described by a set of Langevin equations 
obtained considering the fluctuations around the mean 
values of the operators in the problem and neglecting 
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FIG. 1: Sketch of the system considered. Fields a and b 
interact with a movable mirror. The two cavities are driven 
by input fields with power Pa.b- Input (output) fields are 
indicated as j'" with j = a,b. ^ describes the Brownian 

motion of the mirror at temperature T. 



any resulting non-linear term. This is a well-established 
tool allowing for the exact reconstruction of the quantum 
statistical properties of the system, as far as the fluctu- 
ations of the operators are small compared to the mean 
values .8]. By defining the equilibrium position of the 
mirror = J^j Goj\asj\'^ /ujm, the stationary amplitudes 
of the intracavity fields asj = + AW/^ ^ and 



the effective detunings 



Gojqs, the linearized 



Langevin equations for fluctuations read 



dtSp = -LOmSq - JrnSp 



j=a,b 



dtSxj = -Kjdxj + AjSy-j + {-iy'''SX 



in 
j ' 



dtSyj 



We have introduced the quadrature operators asso- 
ciated with the input noise to the cavities 6Ql^ = 



-y^Kj^gj", (Q — X,Y; q — x,y) and the operator ^ ac- 
counting for the zero-mean Brownian noise. The input 
noise is correlated as {So,in^j{t)5a\^ ^,{t')) — SjkS{t — t') 

with Sdinj — {Sxin.j + iSyinj)/\/2. Moreover, for 
7„ « LOrn we have {i{t)^it') + Ht -t')M- 

The linearity of Eqs. lj2Jl preserves the Gaussian char- 
acter of the system. We can define the vector f-'" = 
{Sxa,Sya,Sxb,6yb,Sq,dp), the kernel K {Gj = V2as^jGoj) 



We aim at studying the entanglement properties of the 
steady state, which is guaranteed to exist if the real parts 
of the eigenvalues of K are negative. For the purposes 
of our work it is sufficient to state that this requirement 
is equivalent to the positivity of two functions, named 
Ci and C2 , the latter of which can be constructed as de- 
scribed in Ref . Ilj . We assume the numbers in the cap- 
tion of Figs.[21for cavity a, which are very close to those 
of recently performed experiments on micromechanical 
systems 6] and, to simplify the calculations, Kb = Ka 
and Gob = Goa (which can be easily relaxed). This al- 
lows us to study Ci,2 as functions of A;, and Pb- We take 
Aa — ujm as this choice corresponds to the maximum 
entanglement between a and the mirror 0] and we con- 
servatively assume that this holds also in presence of b 
(which is a good approximation if P;, <C Pa)- With these 
choices, the behavior of Ci_2 is shown in Figs. HJ Even 
though for Pb < Pa any sign of Af, corresponds to a sta- 
ble regime, we focus on the region associated to A;, < 
as we want to study the interaction of fields a and b with 
the mirror for any value of the back-action induced by b. 

Intracavity entanglement. - At the steady state, 
f ((X)) = fss = limt^oo /q dre^'^hit — t). The stationary 



covariance matrix V, 



PI 



{fss.pfss^q ~t~ fss^qfss^p)/'^ the 



(2) tripartite system can be written as 




(4) 



where Lj accounts for the local properties of subsystem 
j = a,b,m. Cjk describes the correlations between j and 
k. The evaluation of V is performed using the Lyapunov 
equation VK + K V — — N, which is found by noticing 
that, in the Markovian hmit, V = dr(e^^)N(e^^)^. 
The Lyapunov equation is linear in the elements of V, 
which can be easily determined, even though the formal 
solutions are cumbersome. 

We can now study the behavior of the entanglement 
between the elements forming the tripartite system. In 
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matrix 



,{t)ng{t') 



associated with the 



correlation 
ng{t')np{t))/2 = NpgSit - t') 
noise vector {t) = {5X^^,5Y^ 
Here N =Ko]l2 ® Kf,]l2®H where H ~ Diag[0, 7„(2n 1)] 
with n = (e'^'^™ — and I2 the 2x2 identity matrix. 

°K*f(0)-f /p dTeK^n(i-r). 



Eqs. Q are solved as f(t) 




FIG. 2: Ci (panel (a)) and C2 (panel (b)) vs. Ab G 
[-1.5Ka, 1.5Ka] aud Pi, e [0, Pa] for Aa = oJm- The hori- 
zontal plane corresponds to zero and is a help to the eye. We 
used (a;^,w;j,7„,Kj)/27r = (10^ 3.7 x lO", 100, 8.8 x 10^)Hz, 
Ij = 1mm, T = 0.4K, Pa = 50mW and fi = 5ng. 



FIG. 3: Logarithmic negativity S'ff^ (panel (a)) and 
(panel (b)) vs. At £ [— 3ti;„i,0] for the parameters in the 
caption of Fig|5| Each line corresponds to a specific value of 

Pb e [0,Pa]. 



what follows, we characterize and quantify the bipar- 
tite entanglement in each intracavity field-mirror sub- 
system and in the field-field one. Quantitatively, we 
adopt the logarithmic negativity Sj^ 0|, which is an 
entanglement monotone |l3| and can be calculated us- 
ing the symplectic spectrum of the partially transposed 
reduced covariance matrix = PV^feP. Here, Y jk 
is the 4x4 submatrix extracted from V by consider- 
ing the blocks in Eq. Q) relative to subsystems j and 
k only and P = I2 © cr^ (with cr^ the r-Pauli ma- 
trix and r = x,y,z). The symplectic spectum {n±} is 
given by the eigenvalues of |iSVj^.| with S = iay ® ia-y 
Explicitly (n±,)2 ^ [^-^ ± (^-2 _ 4 det V,fc)i/2]/2, 
with xfk = detLj -I- detLfc ± 2detCjfe Q. Entan- 
glement in the state described by Vj^ is found when 
nji^ < 1/2, which translates the criterion for insepa- 
rability of Gaussian states (based on the negativity of 
the partial transposition criterion (NPT) jl^]) in the for- 
malism of the symplectic spectrum. With these tools, 
f^'-=max[0,-ln(2n-,)]. 

We start with the entanglement between field a and 
mirror m. Without field 6, the a — m entanglement 
achieves its maximum 4j. However, the back action in- 
duced by the b field could distort this picture, affecting 
the a — m entanglement. A way to see this is to fix 
the working point for the a cavity to those values cor- 
responding to the maximum of entanglement with the 
mirror. Then, as the effects of the b field are tuned, we 
study the behavior of the a — m entanglement. To this 
task, as done before, we vary Ah and Pt and examine the 
changes in entanglement from negligible to strong back 
action induced by b. The results are shown in Fig. |2| 

As Pb increases, the back-action of field & on m becomes 
more relevant, thus affecting the a — m entanglement at 
moderate values of A^. However, if Af, grows, £^ re- 
vives achieving again values close to its maximum (which 
is larger than 0.3 0]), even more evidently if the range 
of A;, is increased as in Fig. ]E\ Indeed, in a far-detuned 
cavity, less input power enters thus taking back the sys- 
tem to a situation of small back-action. The reduction in 
£^ is caused by a simultaneous raise of S'^, as shown 
in Fig.|31(b). Indeed, the calculation of reveals that 



FIG. 4: Logarithmic negativity £^ vs. Ab€ [~3LJm,0] for the 
parameters in the caption of Fig. |21 Each line corresponds to 
a value of Pb G [0,Pa]- 

entanglement is established between the two subsystems, 
pronounced in the region of moderate A;, where £^ suf- 
fered the effects of 6's back-action. In this case, large 
detunings lower £^"' which, eventually, goes to zero as 

Ab > w„ (Fig.m. 

The complementary behavior of £j^^ 's is an evidence of 
the way the presence of b enables us to infer the features 
of system a — m: the ancillary field b gets entangled with 
a and m, at the expense of the E^""- Given the symmetry 
between a and 6, this same claim holds if we swap o and b. 
Indeed, the most interesting aspect of the entanglement 
dynamics comes from the study of £^ (see Fig. 0J|. As 
soon as is established, £^ becomes non-zero. By 
comparing the results obtained for increasing Pb and A^ 
fixed at the value corresponding to the maximum of £^ , 
we see that £^ disappears, slowly with respect to £^ 
(at Pb ~ Pa, = and ^ 0). 

We may use the entanglement between a and 6 as a 
tool to see signatures of entanglement between a and to, 
as a and b never directly interact and all £^ 7^ is nec- 
essarily due to a mediation by the mirror. That is, in 
this system the mirror acts as a bus for the cross-talking 
of the fields. Any entanglement between a and b is thus 
an indication of a coherent field-field interaction through 
radiation pressure. As the input fields are prepared in 
pure coherent states, £^ must be the result of an effec- 
tive entangling field-field interaction Ifl]. Moreover, by 
finding £^ ^ we can infer entanglement between one of 
the fields and the mirror, even though the converse is not 
true (there are situations where £^ — with £j^" 7^ 0, 
as in Fig. 131(a)). 

Therefore, £^ can be taken as a signature of entangle- 
ment between the cavity a and the mirror to. If initially, 
a and b are in pure states, £^ > strictly indicates en- 
tanglement between a — m. While, in general, one could 
construct models where two systems become entangled 
via the coupling to a system that remains separable with 
respect to the rest |l7j , the interaction studied in our in- 
direct scheme provides strong evidence for mirror-cavity 
entanglement. The case of Fig. El(b) is interesting: for 
Ab/ujm ^ 0.5, - ^aT with achieving its max- 
imum, thus optimizing the overall entanglement distri- 
bution within the system. By means of NPT we have 
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FIG. 5: vs. lAt/tj^l G [0,8] for increasing Pt- Panel (a) 
is for Pi = 0, (b) for Pt = 0.15 Pa and (c) for Pt ~ O.GPa- 
Panel (b) shows a situation where entanglement is found in 
any bipartite system obtained by tracing out one party. 



checked that, in these conditions, genuine tripartite en- 
tanglement is shared between the subsystems. The study 
of an entanglement monogamy inequality in our system 
and lower bounds exploiting a promise to the interaction 
will be the focus of further investigations . 

Extracavity description. - Even though the entangle- 
ment between the intracavity fields and the mirror is 
the object of our investigation, the accessible quanti- 
ties in this system are given by the fields leaking out 
of the cavities. We now show a simple operative strat- 
egy to infer the correlation properties of the intracav- 
ity system. We aim at estimating the field-field covari- 
ance matrix at the output. We assume Kj — k (the 
generalization is straightforward) and define f^{t) = 
fip/y/2K {p = 1, ..,4). The extracavity field quadratures 
put ^ (^g^out^g^out^g^out^g^out^^ related to the in- 
tracavity ones by the input-output relations fp 



The outputs are free fields and 



their dimension is sec~^/^. ft is thus convenient to intro- 
duce dimensionless extracavity quadratures which we use 
to build up the output covariance matrix. One way is to 
define /^^^ = lim^^o, f^{t')df, [u = m.out), 

where tm is the measurement time, i.e., the acquisition 
time chosen for a measurement of the output quadratures 
at the steady state. It is easy to check that /^p's satisfy 
the usual canonical commutation rules. In this way, the 
input-output relations become j^^* — y^2K,t,nfss,p — fj^p, 
where we used /p(oo) = fss.p- 

With this notation, V;^* = {ff^ + fZ^Z)/'^ 
easily evaluated. We find V""* = 2ntjn'Vab + 'V™ , where 



ym 

pq 



i(5 



2^pq- The simplicity of the expression relating 
the extracavity correlations to the analogous intracavity 
quantities suggests an operative way to infer the entan- 
glement behavior of fields a and b. For a fixed working 
point and a value for tm (typically ^ , V°"* is built 
up by homodyne measurements 4, 20]. Then, Vah can be 
reconstructed as (V°"* — V'")/2Kt„i. This prescription is 
just an additional step in the numerical postprocessing of 
the data required for the estimation of the entanglement. 

Conclusions. - We have introduced a scheme to reveal 
entanglement between a cavity field and a movable mirror 
by inducing quantum correlations in the tripartite system 



which includes an ancillary field. Using state of the art 
parameters, we have studied the effects of back-action by 
the ancilla on the entanglement which has to be inferred. 
We found a working point at which entanglement ap- 
pears in any bipartite subsystem. Present work includes 
the formulation of a lower bound to optomechanical en- 
tanglement related to the detected optical one, which will 
allow for the use of field-field correlation as a quantitative 
witness for the mirror-field entanglement. We hope that 
the presented work will pave the way to the experimental 
inference of entanglement involving macroscopic objects. 
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